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Using large scale numerical simulations, we examine the ordering of colloidal particles on square periodic two-dimensional 
muffin-tin substrates consisting of a flat surface with localized pinning sites. We show that when there are four particles per 
pinning site, the particles adopt a hexagonal ordering, while for five particles per pinning site, a square ordering appears. For fill- 
ings between four and five particles per pinning site, we identify a rich variety of distinct ordering regimes, including disordered 
grain boundaries, crystalline stripe structures, superlattice orderings, and disordered patchy arrangements. We characterize the 
different regimes using Voronoi analysis, energy dispersion, and ordering of the domains. We show that many of the boundary 
formation features we observe occur for a wide range of other fillings. Our results demonstrate that grain boundary tailoring can 
be achieved with muffin-tin periodic pinning substrates. 
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1 Introduction 

There is tremendous interest in understanding how to create 
different types of self-assembled colloidal structures. It would 
be particularly useful to identify methods for controlling do- 
main formation and morphology. In two-dimensional systems 
of colloidal particles interacting with a repulsive Yukawa po- 
i tential, the lowest energy configuration is a hexagonal lattice 1 , 
but this ordering can be modified if the particles interact with a 
substrate potential. For random substrates, dislocations appear 
in the colloidal lattice and destroy the hexagonal ordering!2El 
For periodic substrates, the hexagonal ordering may persist, 
but it is also possible for different ordered states as well as 
frustrated partially-ordered states to arise due to the competi- 
tion between the colloidal interactions and the substrate. 

For one-dimensional periodic substrates, triangular and 
smectic type orderings have been observed!^. For two di- 
mensional periodic substrates, the type of ordering depends 
on whether the number of colloidal particles is greater or less 
than the number of potential minima. In the case of egg-carton 
substrates, colloidal molecular crystal states appear whenever 
the number of colloidal particles per potential minima or the 
filling factor / is an integer with / = 2 or higher. Here, the n 
particles in each potential trap form an ft-mer state with an ori- 
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entational degree of freedom, and this system can be mapped 
to various spin modelsE3£Hl Another type of two-dimensional 
substrate that has been experimentally realized is a muffin tin 
where the potential minima or pinning sites are defined as 
localized trapping sites of radius R s surrounded by a flat re- 
gion 15 17 . Experimental realizations of this system show that 
in this case there are two types of colloidal particle species: the 
particles that are directly trapped at the pinning sites, and the 
remaining particles that are localized on the flat part of the po- 
tential due to their interaction with the directly trapped parti- 
cles. These are termed inters titially pinned particles, and they 
are much more mobile than the particles in an egg-carton po- 
tential 15 . Muffin tin-type potentials have been studied for the 
ordering and pinning of vortices in type-II superconductors, 
where different types of vortex crystalline states were shown 
to occur at integer fillings! 1 8 " 2Q i 

At noninteger incommensurate fillings, the system may be- 
come disordered. For colloidal particles in egg-carton poten- 
tials at fillings below the / = 1 matching filling, certain crys- 
talline type structures appear at rational fractional filling ra- 
tios such as / = 1/3, while at other fillings, the system is ei- 
ther disordered or else domain wall structures form 21 . The 
ordering of repulsive particles for fillings less than one has 
been studied for vortices in Josephson junction arrays and in 
periodic pinning, where ordered and quasi-ordered patterns 
ariseGSSIES Such structures also have many strong similari- 
ties to atomic ordering on periodic substrates where the atomic 
coverage is less than a monolayer 25 ! 26 ! The incommensurate 
structures for colloidal particles on muffin type structures at 
higher filling fractions is not known; the additional degrees of 
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Fig. 1 The Voronoi diagrams of the colloidal particle configurations 
on a square pinning array. Filled dots: particle locations; open 
circles: pinning site locations. Polygon colors indicate the 
coordination number of the particle at the center of the polygon: 4 
(dark blue), 5 (light blue), 6 (grey), or 7 (red), (a) At / = 4, each 
pinning site captures one particle and a hexagonal lattice forms, (b) 
At f = 4.5, the particles form a superlattice structure with a 
periodicity of twice the pinning lattice unit cell, (c) At / = 5.0, a 
square lattice forms, (d) A stacked percentage chart of the fraction 
of particles with coordination number of 4 (dark blue), 5 (light 
blue), 6 (grey), and 7 (red) vs filling factor / showing how the 
system evolves from a hexagonal structure at / = 4 to a square 
structure at / = 5. 



2 Simulation and System 

We consider a two-dimensional system with periodic bound- 
ary conditions in x and y containing N c colloidal particles with 
a density p = N c /L 2 , where L is the size of the sample in 
each direction. The particle-particle interaction is given by a 
Yukawa or screened Coulomb potential, and the colloids also 
interact with a square array of N p pinning sites that are each 
of radius R p . The filling fraction is f — N c /N p . The parti- 
cle configurations are obtained by starting from a high tem- 
perature state and performing simulated annealing to reach a 
frozen state. The motion of particle i during the annealing 
arises from integrating the following overdamped equation of 
motion: 

dRi 

dt 



fj^r = -Iw(*y)+FT 

¥J 



(1) 



where 7] is the damping coefficient. The repulsive 
particle-particle Yukawa interaction potential is V(Rij) = 
Eq exp ( — KRij ) / Rij , where R tj = |R; -R/|, R;(y) is the po- 
sition of particle Eq — Z* 2 /(4;r££o<2o)> Z* is the effec- 
tive charge, £ is the solvent dielectric constant, and 1 / K is 
the screening length. The pinning interaction is modeled 
as arising from non-overlapping parabolic traps with Ff = 

F p (R ik /R p )e(R p - R ik )R ik , where R ik = |R f - R*| is the 
distance between particle i and the center of pinning site k lo- 
cated at R^, R/fc = (R; — Rfc) /Rfk, F p is the maximum force of 
the pinning site, and is the Heaviside step function. The ef- 
fects of thermal fluctuations are modeled by the Langevin term 
F T in the form of randomly distributed thermal kicks with the 
properties (F T (t)) = and (F? (t)F(tj)) = 2j]k B 8ij8{t -t f ), 
where k B is the Boltzmann constant. We have tested several 
annealing rates and work with a rate that is slow enough that 
no noticeable difference in the colloidal arrangement appears 
if an even slower rate is used. 



freedom in the interstitial regions could allow for novel order- 
ings or grain boundary structures that cannot form at incom- 
mensurate fillings on egg-carton potentials. 

Here we show that colloidal particles at noninteger fillings 
on muffin tin potentials exhibit a remarkable variety of order- 
ings, including a stripe crystalline regime dominated by five- 
fold and sevenfold coordinated particles. We also find other 
structures including grain boundaries, superlattices, and dis- 
ordered patchy regimes. We specifically focus on fillings be- 
tween four and five particles per trap; however, our results can 
be generalized to many other non-integer fillings with more 
than three particles per trap, and indicate that distinct grain 
boundary morphologies are a generic feature of this system. 



3 Colloidal Particle Configurations 

We concentrate on systems with N p = 400 and N c = 1600 to 
2000, giving a filling fraction 4 < / < 5. At / = 4 each pin- 
ning site captures one particle while the remaining particles 
are located in the interstitial regions. The result is a hexago- 
nal lattice structure in which all of the particles have exactly 
six neighbors. In Fig. [TJa) we plot the Voronoi diagram for 
the / = 4 state indicating the coordination number of each 
particle. Here, a hexagonal lattice forms and all particles are 
sixfold coordinated, while at / = 5.0, Fig.[TJc) shows that the 
particles form a square lattice with fourfold coordination. At 
/ = 4.5 in Fig.[TJb), an exotic isotropic superlattice state ap- 
pears consisting of alternating square arrangements of fivefold 
coordinated particles with a fourfold coordinated particle in 
the center and groups of six sixfold coordinated particles. The 
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(a) f = 4 (b) f = 4.01 5 (c) f = 4.035 (d) f = 4.055 




Fig. 2 The Voronoi diagrams with the same coloring convention as in Fig. 1 for fillings / = (a) 4.0, (b) 4.015, (c) 4.035, (d) 4.055, (e) 4.080, 
(f) 4.085, (g) 4.090, and (h) 4.095. In this disordered domain wall regime, the particles form grain boundaries of paired fivefold and sevenfold 
coordinated (5-7) defects separating different grain orientations. 



unit cell of the structure in Fig. [TJb) is twice the size of the 
pinning lattice unit cell. 

We next address how the system evolves between the hexag- 
onal and square orderings as / is varied from / = 4 to / = 5. 
In Fig. [TJd) we plot the fraction of particles with coordina- 
tion number 4, 5, 6, and 7 for 4.0 < f < 5.0, showing that 
the system begins with hexagonal ordering and has a peak in 
the fraction of fivefold and sevenfold coordinated particles at 
/ = 4.5. For f> 4.65 we observe the rapid growth of the frac- 
tion of fourfold coordinated particles with increasing / until 
the system reaches a square lattice at / = 5.0. We identify four 
regimes as a function of /: a disordered domain wall regime 
for 4.0 < / < 4.1, a stripe crystal regime for 4.1 < f < 4.6, a 
disordered patchy regime for 4.6 < / < 5.0, and the crystalline 
states at / = 4, 4.5, and 5. 

In Fig. [2] we illustrate the domain wall regime for / = 4, 
4.015, 4.035, 4.055, 4.080, 4.085, 4.090, and 4.095. In this 
regime, there are no individual defects; instead, pairs of five- 
fold and sevenfold coordinated (5-7) defects assemble into 
grain boundaries separating two different orientations of the 
hexagonal ordering. The grain boundaries grow in size un- 
til they span the entire system, as shown in Fig. [2jh) for 
f = 4.095. If we repeat the simulations using different ini- 
tial randomizations of the temperature fluctuations, we obtain 
grain boundaries of the same general shape; however, the ex- 
act locations of the grain boundaries can differ. 

In Fig. [3] we show particle configurations in the stripe 
regime for / = 4.1, 4.25, 4.3, 4.35, 4.4, 4.45, 4.5, and 4.55. 



Here the system forms stripe arrangements of 5-7 defects, 
where the number of stripes grows with increasing /, and 
with the isotropic superlattice structure appearing at / = 4.5 
in Fig.[3jg). As / changes, the orientation of the stripes with 
respect to the underlying pinning lattice can change, and the 
thickness of the stripes can vary. Some stripe states consist of 
one-dimensional lines, while in other stripe states, such as at 
/ = 4.1 in Fig.[3ja), there are zig-zag patterns. These patterns 
can even be combined, such as at / = 4.4 in Fig. |3je) where 
there are both zig-zag and one-dimensional lines of stripes. At 
/ = 4.45 and / = 4.55 in Fig.^f) and Fig.^h) there are two 
stripe structures that are interspersed. We find similar stripe 
patterns at these fillings in larger samples. Since the stripe di- 
rections are degenerate, it is possible that domains of different 
stripe orientation could form in the sample, as found for the 
ordering of repulsive particles such as superconducting vor- 
tices on periodic arrays for fillings below / = v 20 i 2l [ In our 
system, since most of the particles are not directly pinned by 
the pinning sites, they can move freely through the intersti- 
tial regions during the annealing process and can more readily 
fall into a global low energy state without different domains. 
In contrast, for / < 1 all the particles are directly pinned by 
the substrate and must hop from one site to another via ther- 
mal activation, producing much stronger kinetic constraints 
and making it more likely that domains of different orienta- 
tion will be quenched into the sample. This result indicates 
that muffin tin potentials permit much larger ordered regions 
than can be obtained with egg-carton potentials. 
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(a) f = 4.1 (b)f = 4.25 (c)f = 4.3 (d)f = 4.35 





Fig. 3 The Voronoi diagrams with the same coloring conventions as in Fig.|l|for fillings / = (a) 4.1, (b) 4.25, (c) 4.3, (d) 4.35, (e) 4.4, (f) 
4.45, (g) 4.5, and (h) 4.55. Here the system forms stripe crystal states containing ordered lines of 5-7 defects, with the number of stripes 
increasing with increasing /. The orientation of the stripes with respect to the underlying pinning lattice varies with filling. At / = 4.5 in 
panel (g), the system forms the isotropic superlattice structure. 






Fig. 4 The Voronoi diagrams with the same coloring conventions as in Fig. 1 for fillings / = (a) 4.65, (b) 4.7, (c) 4.75, (d) 4.8, (e) 4.85, (f) 4.9, 
(g) 4.95, and (h) 5.0. Here the stripe structure is replaced by disordered patchy regions of square domains that grow in size with increasing /. 
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Fig. 5 (a) The fraction of particles in minority grains F m (circles) 
and majority grains Fm (squares) vs /. For / < 4.1 the system is in 
the domain wall regime. When the two curves meet, the system 
enters the stripe regime. For / > 4.6, stripes can no longer be 
identified and the system enters the disordered patch regime, (b) 
Distribution of the stripe orientation angle 6 with respect to the x 
axis vs / in the stripe regime. The stripes are predominantly aligned 
along 6 = 45°, 53°, 59°, and 66.8°, corresponding to symmetry 
directions of the underlying square pinning array. 



Figure]?] shows the evolution out of the stripe state into the 
disordered patch state at / = 4.65, 4.7, 4.75, 4.8, 4.85, 4.9, 
4.95, and 5.0. The stripe state is still present at / = 4.6; how- 
ever, at / = 4.65 in Fig. [4j a )> patches of disorder begin to 
emerge, and at / = 4.70 in Fig. |4jb) some fourfold coordi- 
nated particles appear in patches that increase in size with in- 
creasing / before finally filling the entire system at / = 5.0 in 
Fig.[4jh). The patches do not form regular stripelike patterns 
as seen for f < 4.6 but are more disordered, implying that the 
low energy state in this regime must be much more degenerate 
than in the stripe forming regime. 



4 Analysis 

In order to quantify where the different regimes occur, we 
have developed an algorithm to identify the grain bound- 
aries and their orientation 27 . In the domain wall regime for 

4.0 < f < 4.1, the sixfold coordinated particles in each do- 
main can take one of two degenerate orientations. We define 
the majority grain direction to be the orientation of the largest 
domain in the sample, and the minority grain direction as the 
other degenerate orientation. We then measure the fraction 
of particles Fm that reside in domains with the majority ori- 
entation and the fraction of particles F m in domains with the 
minority orientation. This definition is necessary since the ma- 
jority grain direction can vary from one disorder realization to 
another. Figure [5ja) shows Fm and F m versus /. Just above 
/ = 4, there are small patches of minority grains, with most 
of the particles in one large majority grain. As / increases, 
Fm and F m approach each other until meeting near / = 4.1 
where the stripe regime begins. The stripe regime consists of 
crystalline states that have approximately equal quantities of 
grains of each orientation. The overall fraction of sixfold coor- 
dinated particles drops with increasing / as more stripe walls 
appear. Our grain boundary identification algorithm breaks 
down for / > 4.6 when the disordered patches appear, and 
also at / = 4.5 in the isotropic superlattice. 

In Fig. [5jb) we plot the distribution of the angle 6 be- 
tween the stripe orientation and the x axis for the fillings 

4.1 < / < 4.6 where the stripes can be identified. There are 
four predominant angles, 6 = 45°, 53°, 59°, and 66.8°. These 
angles can be written as 6 — tan -1 (p/q) with integer p and q 
for p/q = 1, 4/3, 5/3, and 7/3, indicating that the stripes are 
aligning with symmetry directions of the underlying square 
pinning array. The other observed values of 6 can also be 
matched to higher order rational ratios of p/q. The alignment 
of particle structures with the symmetry directions of an un- 
derlying pinning lattice has also been observed for colloidal 
particle ordering on quasicrystalline arrays 28 as well as for 
driven colloidal particles moving over periodic and quasicrys- 
talline arrays'SSEI. Not all angles corresponding to possible 
rational values of p/q can be realized due to the fact that the 
domain walls are composed of 5-7 defects and therefore have 
a finite thickness. 

For each filling, we can determine the particle-particle in- 
teraction energy of each particle E{ = Y^jV(Rij) an d then 
construct the distribution P(E) of the energy of all particles 
in the system. In Fig. [6] we plot the resulting P(E) versus 
/. For / < 4.1, P(E) has a maximum just above E = 6.2, 
corresponding to the per-particle energy associated with the 
hexagonal lattice at / = 4.0. As / increases, the peak in P(E) 
remains close to E = 6.2 with some weight in P(E) shifting to 
higher energies due to the increasing fraction of grain bound- 
aries in the sample. The highest energy particles are located 
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Fig. 6 The energy distribution per particle P(E) vs / plotted using a 
color scheme where red indicates high values of P(E) and blue 
indicates low values of P(E). In both the grain boundary regime and 
the disordered patch regime, P(E) is broad, while in the crystalline 
stripe regime, P(E) is much narrower. 



Fig. 7 A stacked percentage chart of the fraction of particles with 
coordination number of 4 (dark blue), 5 (light blue), 6 (grey), and 7 
(red) vs / over the range 1.0 <f < 9.0. Square lattices appear at 
/ = 1.0, 2.0, and 5.0, while hexagonal lattices form at / = 4.0, 6.0, 
8.0, and 9.0. 



along the grain boundaries. Once the system enters the stripe 
regime for 4.1 < f < 4.6, P(E) becomes much narrower due 
to the strong crystalline ordering of the stripe phases, while the 
peak energy increases linearly with increasing E. In the disor- 
dered patch regime, P(E) broadens again, indicating that the 
particles form disordered states with many different particle- 
particle interaction energies. Just below and at / = 5.0, P(E) 
becomes very narrow again as the square lattice forms. 

We can compare the stripe ordering to the pattern forma- 
tion found on periodic egg-carton potentials such as vortex or- 
dering in superconducting wire networks or colloidal particles 
on optical trap arrays for fillings / < I.cP^ES. In these net- 
work systems, the / = 1/3 filling for a square array produces 
a stripelike pattern of particles. This is contrast to our system, 
where the stripes are not individual particles but domains of 
non-sixfold coordinated particles. In the network systems, the 
ordering from < / < 1/2 has a duality with 1/2 < / < 1, 
where the configurations for / = 1/2 — x and / = 1/2 +x are 
the same except that particle locations are replaced by void 
locations. In our system for 4.0 < / < 5.0, a similar duality 
does not occur. Additionally, at the higher fillings, the inter- 
stitial regions in the muffin tin potential allow particles to take 
positions that are inaccessible in egg-carton potentials. 

5 Other fillings 

We have also examined incommensurate fillings from / = 1.0 
to / = 9.0. In Fig. [7] we show the fraction of particles with 



coordination number 4, 5, 6, and 7 versus / over this entire 
range. Square lattices occur at / = 1.0, 2.0, and 5.0, with 
100% of the particles fourfold coordinated at these fillings, 
while hexagonal lattices occur at / = 4.0, 6.0, 8.0, and 9.0, 
with 100% of the particles sixfold coordinated at these fill- 
ings. At / = 7.0 the system forms a patchy nonordered pat- 
tern, while at / = 3.0 the interstitial colloidal particles form a 
dimer crystal arrangement with alternating dimer orientation. 
Figure [7] shows considerable fine structure at noninteger fill- 
ings, and certain types of superlattices form for some fillings 
as indicated by pronounced peaks and dips at / = 3.5, 4.5, 
and 5.5. In general, we find that domain wall formation oc- 
curs at most incommensurate fillings; however, the morphol- 
ogy can exhibit differences depending on the closest integer 
filling. The best examples of stripe crystalline states appear 
between / = 4 and / = 5. In the range / = 2 to / = 3, the sys- 
tem forms only disordered patch states, and no grain bound- 
aries could be identified, while between f = 3 and / = 4, the 
system exhibits grain boundary and stripe regimes but no dis- 
ordered patchy regime. The ordering at / = 3 is not triangu- 
lar or square but is instead a superlattice state. Figure [SJa) 
shows the grain boundary state at / = 3.075, Fig.[8jb) shows 
the stripe state at / = 3.4, and Fig.[8jc) shows that at / = 3.925 
a grain boundary state appears that is similar to that found at 
f = 4.05. For 5 < / < 6, the system exhibits disordered patch 
states of the type illustrated in Fig.[8jd) for / = 5. 1, as well as 
stripe states where the stripes are aligned only along 6 = 0° or 
9 = 90°, as shown in Fig. [SJe) for / = 5.35. For 6 < / < 7, 
the system does not form stripes but is instead disordered, as 
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Fig. 8 The Voronoi diagrams with the same coloring conventions as in Fig.[l]for other fillings, (a) At / = 3.075, grain boundaries appear, (b) 
At / = 3.4, stripes form, (c) At / = 3.925, grain boundaries again form, (d) At / = 5.1 there is a disordered patch state, (e) At / = 5.35 a 
stripe state forms, (f) At / = 6.35 the system is disordered, (g) At / = 7.6, labyrinth type grain boundaries form, (h) At / = 8.8, a particle 
stripe state appears. 



shown in Fig.[SJf) for / = 6.35. For 7 < / < 8, we mostly find 
labyrinth type grain boundary states such as that illustrated in 
Fig.[8jg) for n — 7.6, while for 8 < / < 9, the system mostly 
exhibits isolated defects in a triangular lattice, with some in- 
stances of stripelike domains as shown in Fig.[8jh) for / = 8.8. 

In our study we have only considered square pinning ar- 
rays; however, the nature of the stripe and grain boundary 
states may change considerably for the case of hexagonal pin- 
ning arrays. Additionally, we have only explored the static 
states, but it is likely that the effective friction experienced by 
the particles under an applied drive can strongly affect the or- 
dering of the particles, as has been recently demonstrated for 
domain wall formation in systems of colloidal particles inter- 
acting with periodic pinning arrays near the / = 1 filling^El. 

6 Conclusions 

We have shown that a rich variety of distinct defect patterns 
occur for colloidal particles interacting with a square muffin 
tin potential array. The system forms a hexagonal lattice for 
four colloidal particles per pinning site and a square lattice 
for five colloidal particles per pinning site, and for interme- 
diate fillings we identify a domain wall regime, a crystalline 
stripe regime, and a disordered patch regime. We show how 
these different regimes can be be distinguished with Voronoi 
diagrams, grain orientation measurements, and energy disper- 
sion. These ordered defect patterns arise in the muffin tin 



potential due to the fact that a fraction of the colloidal par- 
ticles are located in the interstitial regions between the pin- 
ning sites and have more mobility than particles on an egg- 
carton potential. We show that the pattern formation of the 
grain boundaries also occurs for a range of other fillings up 
to nine colloidal particles per trap, with the most prominent 
stripe crystalline states forming when there are between four 
and five colloidal particles per trap. Our results may provide a 
new method for tailoring defect structures in colloidal particle 
systems. 
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